library(VennDiagram)
## Loading required package: grid
## Loading required package: futile.logger
library(Seurat)
## Loading required package: ggplot2
## Loading required package: cowplot
##
## Attaching package: 'cowplot'
## The following object is masked from 'package:ggplot2':
##
## ggsave
## Loading required package: Matrix
markers <- read.table('../data/markergenes/markergenes-crossed', header=T)
Marker genes of the intersection of Supra and Peri. Supra has 99 marker genes, Peri 71. They have 9 genes in common.
genes.peri <- markers[markers$cluster == 'Peri.white', 'gene']
genes.supra <- markers[markers$cluster == 'Supra.white', 'gene']
peri.supra <- intersect(genes.peri, genes.supra)
grid.newpage()
draw.pairwise.venn(length(genes.peri), length(genes.supra), length(peri.supra), category=c('Peri', 'Supra'), lty=rep('blank', 2), fill=c('#f67770', '#1fbfc3'), scaled=T, cat.fontfamily = rep("helvetica", 2))
## (polygon[GRID.polygon.11], polygon[GRID.polygon.12], polygon[GRID.polygon.13], polygon[GRID.polygon.14], text[GRID.text.15], text[GRID.text.16], text[GRID.text.17], text[GRID.text.18], text[GRID.text.19])
Marker genes of the intersection of Subq and Visce. Subq has 65 marker genes, Visce 22. They have only 2 genes in common.
genes.visce <- markers[markers$cluster == 'Visce.brown', 'gene']
genes.subq <- markers[markers$cluster == 'Subq.brown', 'gene']
visce.subq <- intersect(genes.visce, genes.subq)
grid.newpage()
draw.pairwise.venn(length(genes.visce), length(genes.subq), length(visce.subq), category=c('Visce', 'Subq'), lty=rep('blank', 2), fill=c('#c680fc', '#7dac1f'), scaled=T, cat.fontfamily = rep("helvetica", 2))
## (polygon[GRID.polygon.20], polygon[GRID.polygon.21], polygon[GRID.polygon.22], polygon[GRID.polygon.23], text[GRID.text.24], text[GRID.text.25], text[GRID.text.26], lines[GRID.lines.27], text[GRID.text.28], text[GRID.text.29])
all10x <- readRDS('../data/10x')
print(peri.supra)
## [1] "TM4SF1" "MTRNR2L8" "LY6K" "MGST3" "FHL2" "COL6A2"
## [7] "TAGLN" "HSPA8" "PENK"
markers[which(markers$gene %in% peri.supra),]
VlnPlot(all10x, features.plot=toupper(c('TM4SF1', 'MTRNR2L8', 'LY6K', 'MGST3', 'FHL2', 'COL6A2', 'TAGLN', 'HSPA8', 'PENK')), group.by='sample_name', point.size.use=-1, nCol=2, x.lab.rot=T, size.x.use=10)
print(visce.subq)
## [1] "GREM1" "HOXB7"
print(markers[which(markers$gene %in% visce.subq),])
## cluster p_val avg_logFC pct.1 pct.2 p_val_adj gene
## 190 Visce.brown 0 0.3825937 0.771 0.691 0 GREM1
## 205 Visce.brown 0 0.2671629 0.540 0.324 0 HOXB7
## 220 Subq.brown 0 0.5615335 0.888 0.691 0 GREM1
## 252 Subq.brown 0 0.2726722 0.601 0.324 0 HOXB7
VlnPlot(all10x, features.plot=toupper(c('GREM1', 'HOXB7')), group.by='sample_name', point.size.use=-1, nCol=2, x.lab.rot=T, size.x.use=10)
print(markers[markers$cluster == 'Peri.white',][1:10,])
## cluster p_val avg_logFC pct.1 pct.2 p_val_adj gene
## 1 Peri.white 0 0.8670990 0.928 0.407 0 TM4SF1
## 2 Peri.white 0 0.8098424 0.963 0.681 0 SRGN
## 3 Peri.white 0 0.6577277 0.564 0.161 0 NEFM
## 4 Peri.white 0 0.6504400 0.830 0.468 0 MEST
## 5 Peri.white 0 0.6249721 0.932 0.628 0 RGS4
## 6 Peri.white 0 0.6177285 0.961 0.817 0 PTX3
## 7 Peri.white 0 0.5764050 0.698 0.354 0 MTRNR2L1
## 8 Peri.white 0 0.5387147 0.967 0.802 0 GLIPR1
## 9 Peri.white 0 0.5092244 0.994 0.964 0 CRYAB
## 10 Peri.white 0 0.5058522 0.917 0.656 0 CTSC
VlnPlot(all10x, features.plot=as.vector(markers[markers$cluster == 'Peri.white',][1:10,'gene']), group.by='sample_name', point.size.use=-1, nCol=2, x.lab.rot=T, size.x.use=10)
print(markers[markers$cluster == 'Supra.white',][1:10,])
## cluster p_val avg_logFC pct.1 pct.2 p_val_adj gene
## 81 Supra.white 0 1.7650025 0.786 0.226 0 IGFBP5
## 82 Supra.white 0 0.9351586 0.977 0.832 0 IGFBP3
## 83 Supra.white 0 0.7779768 0.978 0.800 0 AKAP12
## 84 Supra.white 0 0.6136319 0.804 0.475 0 TNFRSF11B
## 85 Supra.white 0 0.5950695 0.955 0.719 0 DCN
## 86 Supra.white 0 0.5918296 0.884 0.539 0 IFI27
## 87 Supra.white 0 0.5818591 0.808 0.448 0 EFEMP1
## 88 Supra.white 0 0.5744373 0.501 0.176 0 PENK
## 89 Supra.white 0 0.5733726 0.995 0.892 0 MFAP5
## 90 Supra.white 0 0.5585209 0.905 0.656 0 CYP1B1
VlnPlot(all10x, features.plot=as.vector(markers[markers$cluster == 'Supra.white',][1:10,'gene']), group.by='sample_name', point.size.use=-1, nCol=2, x.lab.rot=T, size.x.use=10)
print(markers[markers$cluster == 'Visce.brown',][1:10,])
## cluster p_val avg_logFC pct.1 pct.2 p_val_adj gene
## 189 Visce.brown 0 0.4954102 1.000 0.999 0 FTH1
## 190 Visce.brown 0 0.3825937 0.771 0.691 0 GREM1
## 191 Visce.brown 0 0.3562843 1.000 1.000 0 FTL
## 192 Visce.brown 0 0.3431171 0.588 0.457 0 LUM
## 193 Visce.brown 0 0.3366891 0.999 0.996 0 GUK1
## 194 Visce.brown 0 0.3158248 0.595 0.430 0 RP11-173B14.5
## 195 Visce.brown 0 0.3156858 1.000 0.998 0 RPS2
## 196 Visce.brown 0 0.3098157 0.930 0.871 0 PLAC9
## 197 Visce.brown 0 0.3089736 0.546 0.365 0 HOXA5
## 198 Visce.brown 0 0.3059777 0.992 0.985 0 IGFBP4
Interestingly, some genes seem to be found because Supra4 and Subq4 are different.
VlnPlot(all10x, features.plot=as.vector(markers[markers$cluster == 'Visce.brown',][1:10,'gene']), group.by='sample_name', point.size.use=-1, nCol=2, x.lab.rot=T, size.x.use=10)
print(markers[markers$cluster == 'Subq.brown',][1:10,])
## cluster p_val avg_logFC pct.1 pct.2 p_val_adj gene
## 213 Subq.brown 0 1.1149278 0.897 0.722 0 THBS1
## 214 Subq.brown 0 0.9239007 0.974 0.948 0 TIMP3
## 215 Subq.brown 0 0.8931719 1.000 1.000 0 RPS29
## 216 Subq.brown 0 0.8099464 0.402 0.094 0 BCYRN1
## 217 Subq.brown 0 0.7024385 0.831 0.729 0 DKK1
## 218 Subq.brown 0 0.6638062 0.975 0.980 0 RPL36A
## 219 Subq.brown 0 0.6095757 0.973 0.961 0 ATP5I
## 220 Subq.brown 0 0.5615335 0.888 0.691 0 GREM1
## 221 Subq.brown 0 0.5231598 0.654 0.644 0 AC009501.4
## 222 Subq.brown 0 0.5113614 0.997 0.996 0 RPS10
Most of the markers for Subq seem to be found because Subq_3 stands out (RPS29, BCYRN1, RPL36A, ATP5I, AC009501.4, RPS10). Mostly ribosomal genes.
BCYRN1 = brain cytoplasmic RNA
VlnPlot(all10x, features.plot=as.vector(markers[markers$cluster == 'Subq.brown',][1:10,'gene']), group.by='sample_name', point.size.use=-1, nCol=2, x.lab.rot=T, size.x.use=10)